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Abstract 

We report a measurement of the cross section of single top-quark production in the f-channel 
using 1.04 fb"' of pp collision data at ^ 7 TeV recorded with the ATLAS detector at the LHC. 
Selected events feature one electron or muon, missing transverse momentum, and two or three jets, 
exactly one of them identified as originating from a b quark. The cross section is measured by fitting 
the distribution of a multivariate discriminant constructed with a neural network, yielding cr, = 83 ± 
4(stat.);^j"(syst.)pb, which is in good agreement with the prediction of the Standard Model. Using the 
ratio of the measured to the theoretically predicted cross section and assuming that the top-quark- 
related CKM matrix elements obey the relation \V,b\ » the coupling strength at the W-t-b 
vertex is determined to be \V,b\ ^ 1.13:^q [3. assumed that \V,b\ < l a lower limit of \Vtb\ > 0.75 is 
obtained at the 95% confidence level. 



Measurement of the t-channel single top-quark production cross section in pp 
collisions dX ^[s = 1 TeV with the ATLAS detector 

The ATLAS Collaboration 



Abstract 

We report a measurement of the cross section of single top-quark production in the f-channel using 1.04 fb"' of pp 
collision data ai yfs - 1 TeV recorded with the ATLAS detector at the LHC. Selected events feature one electron 
or muon, missing transverse momentum, and two or three jets, exactly one of them identified as originating from 
a b quark. The cross section is measured by fitting the distribution of a multivariate discriminant constructed with 
a neural network, yielding cr, = 83 ± 4 (stat.) ;^ jg (syst.) pb, which is in good agreement with the prediction of the 
Standard Model. Using the ratio of the measured to the theoretically predicted cross section and assuming that the 
top-quark-related CKM matrix elements obey the relation \Vth\ » \Vis\, \Vtd\, the coupling strength at the W-t-b vertex 
is determined to be \V,i,\ - LlS^^ jj. If it is assumed that \Vth\ < 1 a lower limit of \Vtb\ > 0.75 is obtained at the 95% 
confidence level. 



1. Introduction 

At hadron colliders top quarks are predomi- 
nantly produced in pairs (top-antitop) via the flavour- 
conserving strong interaction. Alternative production 
modes proceed via the weak interaction involving a W- 
t-b vertex, leading to a single top-quark intermediate 
state. Three subprocesses contribute to single top-quark 
production: the exchange of a virtual W boson in the 
f-channel, or in the i-channel, and the associated pro- 
duction of a top quark and an on-shell W boson. The 
process with the highest cross section at the Tevatron 
and at the LHC is the f-channel mode q + b ^ q' + t. 

In 2009, single top-quark production was observed 
by flie CDF II and D0 Collaborations based on 
analyses counting the f-channel and .s-channel processes 
as signal. The observation of the f-channel production 
mode has also been recently reported by D0 fl, while 
the CMS Collaboration has published evidence of this 
process at the LHC [4]. 

The single top-quark final state provides a direct 
probe of the W-t-b coupling and is sensitive to many 
models of new physics [5i\. The measurement of the 
production cross section constrains the absolute value of 
the quark-mixing matrix element V,b 7] without as- 
sumptions about the number of quark generations (see 
Ref. ^ for a recent measurement from the D0 Col- 
laboration). Alternatively, it allows the /j-quark parton 
distribution function (PDF) to be measured. 

At the LHC, colliding protons at ^/s - 1 TeV, the 
sum of f and f cross sections is predicted to be: cr, - 
64.6^^ Q pb I9J for the leading f-channel process, ctw, - 
15.7 ± 1.1 pb iH for Wt associated production, and 
CTj = 4.6 + 0.2 pb ifllll for the i-channel. The analyses 



presented in this Letter consider only the f-channel pro- 
cess as signal, while the other two single top-quark pro- 
cesses are treated as backgrounds, assuming the Stan- 
dard Model (SM) theoretical cross sections for these 
processes. 

The W boson from the top-quark decay is recon- 
structed in its leptonic decay modes ev, yuv or rv, where 
the T decays leptonically. Thus, selected events con- 
tain one charged lepton candidate, e or yu; two or three 
hadronic high-px jets; and missing transverse momen- 
tum £'™'^\ Two jets are expected from the leading-order 
(LO) process, while a third jet may arise from higher- 
order processes. Exactly one of the jets is required to be 
identified as originating from a ^j-quark. 

The measurement of cr, is based on a fit to a multi- 
variate discriminant constructed with a neural network 
(NN) to separate signal from background and the result 
is cross-checked using a cut-based method, which ad- 
ditionally provides a breakdown for the f and f cross 
sections. 

2. Data and simulated event samples 

The analyses described in this Letter use proton- 
proton LHC collision data at a centre-of-mass energy 
of 7 TeV collected with the ATLAS detector [12] be- 
tween March and June 2011. The selected events were 
recorded based on single electron and muon triggers. 
Stringent detector and data quality requirements are ap- 
plied, resulting in a data set corresponding to an inte- 
grated luminosity of 1 .04 ± 0.04 fb ' lfllll4l. 

Samples of simulated events for all three single top- 
quark processes are produced with the AcerMC pro- 
gram (version 3.7) using MRST 2007LO* paiton 
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distribution functions (PDFs) III6II . The computation of 
the f-channel single top-quark process in AcerMC in- 
corporates the ^+/7 — > ^'-Hf and^+^ — » ^'+f-i-l7 subpro- 
cesses and features an automated procedure to remove 
the overlap in phase space between them |17]. Sam- 
ples of the top-quark pair (tt) process are generated us- 
ing MC@NLO (version 3.41) [18], with the CTEQ6.6 
set of PDFs 1 19]. The top-quark mass is assumed to 
be 172.5 GeV. Generator default values of 0.999105 
and 0.999152 ai-e used for \ V,b\ to produce the AcerMC 
andMC@NLO samples, respectively. At higher orders 
in perturbation theory, interference effects between the 
single-top Wt channel and tt processes occur, but are 
found to be small [|20 1 and can therefore be safely ne- 
glected. The ALPGEN leading-order generator (version 
2.13) 12?] and the CTEQ6L1 set of PDFs |19] are used 
to generate W-njets, Wbb, Wcc, Wc and Z-i-jets events 
with up to five additional partons. To remove overlaps 
between the n and n + 1 parton samples the "MLM" 
matching scheme 112 ill is used. The double count- 
ing between the inclusive W + n parton samples and 
samples with associated heavy-quark pair-production 
is removed utilising an overlap removal based on a 
AR - -^/(At/)^ + (A0)2 matching. The diboson pro- 
cesses WW, WZ and ZZ are generated using HERWIG 
(version 6.5.20) [.221 . For all single top-quark sam- 
ples the hadronisation is performed by PYTHIA (ver- 
sion 6.4.25) |23]; in all other cases HERWIG in con- 
nection with the JIMMY |24] underlying event model 
(version 4.31) is used. After the event generation, all 
samples are passed through the full simulation of the 
ATLAS detector |25] based on GEANT4 |26] and are 
then reconstructed using the same procedure as for col- 
lision data. The simulation includes the effect of mul- 
tiple pp collisions per bunch crossing and is weighted 
to the same distribution as observed in the data with 5.6 
interactions per bunch crossing on average. 

3. Object definition and event selection 

Electron candidates are reconstructed offline using a 
cluster-based algorithm and are required to have Ej > 
25 GeV and |7/ci| < 2.47, where rj^i denotes the pseudo- 
rapidity of the calorimeter cluster. Clusters in the transi- 
tion regions between the calorimeter barrel and endcaps. 



ATLAS uses a right-handed coordinate system with its origin at 
the nominal interaction point in the centre of the detector and the z- 
axis along the beam direction. The A-axis points towards the centre 
of the LHC ring, the v-axis points upwards. The pseudorapidity is 
defined in terms of the polar angle 6 slst] = - lntan(6/2). 



corresponding to 1.37 < |77cil < 1.52, are ignored. High- 
quahty electron candidates are selected using a set of 
cuts 1I27II which include stringent requirements on the 
matching between the track and the calorimeter clus- 
ter. Electrons must also be isolated: the sum of the 
calorimeter transverse energy within a cone of radius 
AR = 0.3 (excluding the cells associated with the elec- 
tron) must be less than 15% of the electron E-y, and the 
Pt of all tracks within the same cone radius around the 
electron direction, again excluding the track associated 
to the electron, must be less than 10% of the electron 
Ej. 

Muon candidates are reconstructed by combining 
track segments found in the inner detector and in the 
muon spectrometer We only consider muon candidates 
that have pj > 25 GeV and I77I < 2.5. Selected muons 
must additionally satisfy a series of cuts on the num- 
ber of track hits present in the various tracking sub- 
detectors (2§]. Muon candidates are required to be iso- 
lated using the equivalent criteria as applied to electron 
candidates. 

Jets are reconstructed using the anti-A;, algorithm l! 29ll 
with a radius parameter of 0.4, using clusters of adja- 
cent calorimeter cells ll30ll as inputs to the jet cluster- 
ing. The response of the calorimeter is corrected by 
Pt;- and ?7-dependent factors 131] , which are applied 
to each jet to provide an average energy scale correc- 
tion. Jets overlapping with selected electron candidates 
within AR < 0.2 are removed, as in these cases the jet 
and the electron are very likely to correspond to the 
same physics object. Only jets having px > 25 GeV 
and I77I < 4.5 are considered. Jets originating from bot- 
tom quarks are tagged in the region I77I < 2.5 by recon- 
structing secondary and tertiary vertices from the tracks 
associated with each jet and combining lifetime-related 
information with a NN ll32fl . A threshold is applied to 
the /7-tagging algorithm output corresponding to a b- 
tagging efliciency of about 57% and a light-quark jet 
rejection factor (the reciprocal of the efficiency to b-tag 
light quarks) of about 520 for jets in events. The E™"** 
is calculated using clusters of adjacent calorimeter cells 
and corrected for the presence of electrons, muons, and 
jets Q. 

Events are selected if they contain at least one good 
primary vertex candidate ||34|| with a minimum of five 
associated tracks each with pj > 400 MeV. Events con- 
taining jets failing quality criteria [35] are rejected. 

The event selection requires exactly one charged lep- 
ton, e or jj, exactly two or three jets, and £'™^'* > 
25 GeV. A trigger matching requirement is applied 
where the lepton must lie within AR < 0.15 of its 
trigger-level object. Since the multijet background is 
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difficult to model precisely, its contribution is addition- 
ally reduced through a requirement on the transverse 
mass of the lepton-^™"' systenjl: nijiW) > (60 GeV - 

The following samples are defined for this analysis: 
a "/j-tagged sample" with two jets or three jets, exactly 
one of which is Zj-tagged, and a "pretag sample" with 
two or three jets, without making any /7-tagging require- 
ment. We also use a sample containing exactly one b- 
tagged jet to estimate the W-njets flavour composition. 

4. Background estimation 

A large background to the single top-quark final state 
comes from QCD-produced multijet events in which ei- 
ther one of the jets is misidentified as an isolated lep- 
ton or a non-prompt lepton (for example from a /7-quark 
semileptonic decay) appears isolated. Other significant 
backgrounds originate from W^-boson production in as- 
sociation with jets and tt production. Smaller back- 
grounds come from Z-i-jets, Wf-channel and s-channel 
single top-quark production, and diboson production. 
These smaller backgrounds and the tt background are 
normalised to their theoretical predictions. For Z-i-jets 
background the inclusive cross sections are calculated to 
next-to-next-to-leading order (NNLO) with FEWZ (ver- 
sion from March 15, 2009) 1,37.1 . The diboson cross sec- 
tions are normalised to next-to-leading order theoretical 
calculations |38]. The cross section is normalised to 
the approximate NNLO-predicted value obtained using 
HATHOR (version 1.2) Q. 

The multijet background normalisation is obtained 
using a binned maximum-likelihood fit to the E™*** dis- 
tribution in the data, before the application of the 
£miss ^^j.^ using a data-derived template for the multi- 
jet background and templates from Monte-Carlo simu- 
lation for all other processes (top quark, W/Z-njets, di- 
bosons). The multijet template is created using collision 
events that are triggered by a single low-pj jet. Sev- 
eral prescaled trigger streams with different thresh- 
olds are used for that purpose. In the offline selection 
of these events the electron requirement is replaced by a 
jet requirement (jet-electron model). This jet must have 
Pt > 25 GeV, the same acceptance in I77I as the signal 
electron, and (80-95)% of the jet energy deposited in 
the electromagnetic section of the calorimeter. The last 
requirement ensures the orthogonality of the jet-electron 



^ Defined as 

where C denotes the lepton. 



data set to the sample of events with electron candidates 
which feature an electromagnetic energy fraction larger 
than 95%. The jet must also contain at least four tracks, 
thus reducing the contribution from converted photons. 
When selecting the jet-electron sample, events contain- 
ing electron or muon candidates in addition to the jet- 
electron are vetoed. The same model is also used in 
the muon channel. A systematic uncertainty of 50% on 
the multijet background rates was estimated by study- 
ing the impact of pile-up events on the fit results and 
by performing likelihood fits on the m-YiW) distribution. 
The jet-electron model is also used to model the shape 
of kinematic distributions of the multijet background. 

The kinematic distributions of the W+jets back- 
ground, which comprises contributions from W-nheavy 
flavour jets (Wbb+jets, Wcc-i-jets and WcH-jets) and 
W-nlight jets, are taken from samples of simulated 
events, while the normalisation of the flavour compo- 
sition is derived from data. The NN analysis simultane- 
ously determines the normalisation of the W-nlight jets 
and W+ heavy flavour processes when fitting the NN 
discriminant distribution to measure the f-channel sin- 
gle top-quark rate. The cut-based analysis derives nor- 
malisation factors for the W-i-jets processes using the 
event yields in the 1-jet /j-tagged, 2-jet pretag, and 2-jet 
/7-tagged sample, excluding events selected by the cuts 
defined in Section |5] Since the 2-jet /7-tagged sample 
includes some f-channel signal events, despite requiring 
that the events fail the selection of the cut-based anal- 
ysis, an uncertainty of 100% on the expected f-channel 
single top-quark rate is assumed in this normalisation 
procedure. Both estimates of the W+jets backgrounds, 
the one of the NN analysis and the one of the cut-based 
analysis, are in very good agreement with each other 

The predicted and observed event yields, after the ap- 
plication of the selections described in Section [3] are 
given in Table [1] separately for the electron and muon 
channels, in the 2-jet and 3-jet /7-tagged samples. For 
the purpose of this table and the histograms of kine- 
matic distributions the contributions of the tt, Wt and s- 
channel processes have been grouped together into one 
category called "ff, other top". In the subsequent analy- 
sis the electron and muon channels are combined. 

5. Signal and background discrimination 

To separate f-channel single top-quark signal events 
from background several kinematic variables are com- 
bined into one discriminant by employing a NN, that 
also exploits correlations between the variables. The re- 
sult of the NN analysis is corroborated by a cut-based 
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Table 1: Predicted and observed event yields, after selection, in the electron and muon 2-jet and 3-jet fo-tagged samples. The multijet event yields 
ai'e determined with data-driven technique. Contributions from W+jets events are normalised to observed data in control regions as used in the 
cut-based analysis. The uncertainties on the multijet and the W+jets yields are also estimated from data (see Section|6). All other backgrounds and 
the f-channel signal expectation are normalised to theoretical cross sections. Uncertainties on these predictions are only reflecting the uncertainties 
on the theoretical cross section prediction and do not include experimental uncertainties (such as the jet energy scale uncertainty, etc.). 





Electron 


Muon 




2-jet 


3-jet 


2-jet 


3-jet 


single-top f-channel 


447 ± 11 


297 ±7 


492 ± 12 


323 + 8 


tt, other top 


785 ± 52 


1700+120 


801+53 


1740+130 


W-i-light jets 


350+ 100 


128 + 56 


510 + 150 


209 + 91 


W-Hheavy flavour jets 


2600 + 740 


1100 + 400 


3130 + 880 


1270 + 480 


Z-i-jets, diboson 


158 + 63 


96+44 


166+61 


80 + 31 


Multijet 


710+350 


580 + 290 


440+220 


270 + 140 


Total expected 


5050 ± 830 


3900 + 520 


5530 + 930 


3900 + 520 


Data 


5021 


3592 


5592 


3915 



analysis that applies additional criteria to the basic se- 
lection described in Section|3] 



Neural network 
roBayes S 



based discriminant. The Neu- 
41|] tool (version 3.3) is used for 



preprocessing the input variables and for the training of 
the NN. A large number of input variables is studied, 
but only the highest-ranking variables are chosen for 
the training of the NN. The ranking of variables is 
automatically determined as part of the preprocessing 
step and is independent of the training procedure. The 
total correlation k^°^^ of a set of variables to the target 
function, that assumes the value 1 for signal and for 
background events, is computed as a measure of the 
discrimination power of these variables. In an iterative 
procedure, the variables are sorted according to the 
loss in /f™"' that is induced due to their removal from 
the set. Considering the number of simulated events 
used to determine the ranking, one can compute the 
significance of the information loss caused by the 
removal of a certain variable. For the training of the 
NN we use only variables that contribute with more 
than 20 (10) standard deviations to a-J"'"' in the 2-jet 
(3-jet) data set. This choice is a compromise between 
the achievable discrimination power, that increases 
with the number of variables, and the practical aim of 
keeping the number of variables at a manageable level. 

As a result of this optimization procedure 12 kine- 
matic variables are identified that serve as inputs to the 
NN in the 2-jet data set. The most discriminating vari- 
able is the invariant mass of the /^-tagged jet, the charged 
lepton, and the neutrino, m(£vb), which is an estimator 
for the top-quark mass for signal events. In this calcu- 



lation the transverse momentum of the neutrino is given 
by the x- and y-components of the £'™'*^ vector, while 
the unmeasured z-component of the neutrino momen- 
tum, Pz(v), is inferred by imposing a W-boson mass 
constraint on the lepton-neutrino system. Since the con- 
straint leads to a quadratic equation for Pziv), a two- 
fold ambiguity arises. In the case of two real solutions, 
the one with the smaller \p~\ is chosen. If the solutions 
are complex, those are avoided by a kinematic fit that 
rescales the neutrino and py such that the imaginary 
radical vanishes, but keeps the transverse components 
of the neutrino as close as possible to the £'™'*^ The 
second and third most discriminating variables are the 
absolute value of the pseudorapidity of the highest pj 
untagged jet \ri{ju)\ and the transverse energy of the un- 
tagged jet Ej (ju)- Other variables used by the NN in the 
2-jet data set are: the absolute value of A?; between the 
^-tagged jet and the reconstructed W boson lArjib, W)\; 
the absolute value of At] between the /j-tagged jet and 
the highest pj untagged jet {Arjib, the transverse 
momentum of the charged lepton piit); the scalar sum 
of the transverse momenta of the lepton, jets, and 
Hi{(,]tts,, £'™'*^); m-Y(W); the pseudorapidity of the lep- 
ton ri(€); the invariant mass of the ^-tagged jet m{b); 
E™^^; and the invariant mass of the untagged jet and the 
Z7-taggedjet m(jy,b). 

For events with three jets 18 variables are used, the 
most discriminating ones being the invariant mass of the 
two leading jets, m{jij2), m{(vb), and the absolute value 
of the difference in the pseudorapidity of the leading and 
lowest Pt jet, 1^^(71,73)1- Figure [1] shows distributions 
of some of the most discriminating variables in the b- 
tagged 2-jet or 3-jet samples, used in both the NN anal- 
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Figure 1: Discriminating variables in the fo-tagged sample for 2-jet events and 3-jet events. IVIultijet event yields are determined with data-driven 
techniques. Contributions from W+jets events are derived from simulation and normalised to data in control regions with the method employed in 
the cut-based analysis. All other backgrounds and the f-channel signal expectation are nomialised to theoretical cross sections. The last histogram 
bin includes overflows. The figures shown are for 2-jet or 3-jet events, respectively: (a), (e) the invariant mass of the fo-tagged jet, the charged 
lepton, and the neutrino; (b), (f) the scalar sum of the transverse momenta of the lepton, the jets, and E™^"; (c), (g) the absolute value of the 
pseudorapidity of the highest px untagged jet. For 2-jet events Figure (d) shows the absolute value of A;; between the fo-tagged jet and the highest 
Pt untagged jet; and for 3-jet events Figure (h) displays the invariant mass of the three selected jets. 
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Figure 2: Additional discriminating variables used in the NN analysis for 2-jet and 3-jet events. The rate of multijet events is normahzed to the 
estimate obtained from the fit to the E'™^^ distributions. All other component distributions are normalised to the result of the maximum-likelihood 
fit of the NN output. The last histogram bin includes overflows. The figures shown are, for 2-jet events: (a) invariant mass of the highest pj 
untagged jet and the fc-tagged jet; (b) transverse mass of the lepton-i?^"'''' system; and for 3-jet events: (c) invariant mass of the two leading jets; (d) 
absolute value of A;; between the leading jet and the lowest pj jet. 



ysis and the cut-based analysis. The variable m(jij2j3) 
denotes the invariant mass of all selected jets in the 3-jet 
data set. Distributions of additional variables used only 
in the NN approach are shown in Figure |2] 

The agreement between the background model and 
collision data is tested in the large pretag sample for 
each input variable used in the analysis, for various ad- 
ditional control variables, and the NN output distribu- 
tions, which are shown in Figs. |3(a)| and |(b)| In this 
control sample, where the /7-tagging algorithm has not 
yet been applied, the ^-tagged jet is substituted by the 
most central jet, with the requirement that it is within 
I77I < 2.5. Good agreement is found overall, except for 
the It/I distribution of the jet with the highest I77I in the 
pretag data set for which an additional systematic mod- 
elling uncertainty is taken into account (see Section|6]l. 

The NeuroBayes tool combines a three-layer feed- 
forward NN with a complex preprocessing of the input 
variables. By transforming the variables in the prepro- 
cessing step the influence of outliers is largely reduced 
and statistical fluctuations are damped. NeuroBayes ap- 



plies Bayesian regularisation techniques for the train- 
ing process to damp statistical fluctuations in the train- 
ing sample and to avoid overtraining. A certain fraction 
of simulated events (20%) is not included in the train- 
ing sample and is used as an independent test sample to 
check that there is no overtraining. The ratio of signal 
to background events in the training is chosen to be 1:1, 
while the different background processes are weighted 
according to the number of expected events. 

To extract the signal content of the selected sam- 
ple a maximum-likelihood fit is performed to the com- 
plete NN output distributions in the 2-jet and 3-jet data 
sets (see Section |7]i. Fitting all bins of the distribution 
has the advantage of making maximal use of the signal 
events remaining after the event selection, and also al- 
lows the background rates to be constrained by the data. 
The sensitivity to the background rates is given by the 
background dominated region close to zero. The ob- 
served NN output distributions scaled to the fit result 
are shown in Figs.| 3Ic)] and[ 3Id)| for /7-tagged events with 
two or three jets, respectively. 
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Figure 3: (a) and (b): Neural network output distribution for the pretag sample, including the JES uncertainty on the prediction (hatched region). 
The multijet component is normalized to the estimate obtained from the fit to the E'™^^ distributions. All other components are normalised such that 
the total number of expected events in the pretag sample is equal to the observed number of events. The ratio beween the data and the total predicted 
distributions is also shown, (c) and (d): NN output distribution for the 2-jet and 3-jet fo-tagged samples, respectively. All component distributions 
are normalised to the result of the maximum-likelihood fit, except for the component of multijet events that is normalized to the estimate obtained 
from the fit to the £5"*'* distributions. 



Cut-based selection. In the cut-based analysis addi- 
tional selections are applied to a subset of five variables 
used by the NN analysis: |;70u)| > 2, i/T(^ Jets, E^''') > 
210 GeV, and 150 GeV < m(ivb) < 190 GeV. The 2-jet 
selection requires \AT]{b,ju)\ > 1, while the 3-jet selec- 
tion requires that mO'iy^ys) is higher than 450 GeV, to 
further reduce the large tt contribution in this channel. 
The selection cuts were chosen in order to increase the 
expected significance of the f-channel single top-quark 
signal, taking into account systematic uncertainties on 
the background estimate 04211 . 



The positive and negative lepton-charge samples 
are considered separately, as more single-top quark f- 



channel events are expected in the e^/ju^ samples than 
in the e^ jji' samples due to the dominance of valence 
u quarks in the proton over d quarks. The 2-jet and 3-jet 
data sets are also considered separately. 

The signal and background event counts for the pos- 
itive and negative lepton-charge samples are given in 
Table |2] The observed event yields are consistent with 
the SM expectation in each channel. Figure j ^fa)] shows 
the distribution of the lepton charge for events with two 
or three jets after the application of all cut-based selec- 
tions. Figures | ^b)| and| 3Jc)| display the variable mitvb) 
for 2-jet and 3-jet events respectively after applying all 
selections except for the cut on m{{vb). In these fig- 
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Table 2: Event yield for the 2-jet and 3-jet fc-tagged positive and negative lepton-charge channels after the cut-based selection. The multijet and 
W+jets backgrounds are normalised to observed data in control regions, all other samples are normalised to theory cross sections. Uncertainties 
shown include all sources of systematic errors, summed quadratically and without taking into account possible anticorrelations between systematic 
sources and between processes. 





Cut-based 2-jet 


Cut-based 3-jet 




Lepton -H 


Lepton - 


Lepton + 


Lepton - 


single-top f-channel 


85 + 29 


39 ±13 


33.6 ±7.0 


14.6 ±6.2 


tt, other top 


14.0 ±6.4 


12.8+4.2 


10.5 ±4.2 


10.7±7.9 


W-i- light jets 


3.3 + 1.9 


2.0 + 1.2 


0.8 ±1.3 


0.3 ±0.3 


W-i-heavy flavour jets 


39 + 11 


27.1 ±7.5 


8.7 ±6.0 


3.4±3.1 


Z+jets, diboson 


1.1+0.8 


1.0±0.8 


0.3 ±0.2 


0.2 ±0.3 


Multijet 


0.2 + 0.2 


0.3 ±0.3 


1.5±1.1 


3.1 ±2.0 


Total expected 


143 + 31 


83 ±16 


56 ±10 


32 ±11 


S/B 


1.5 


0.9 


1.6 


1.0 


Data 


193 


101 


53 


39 



ures, the f-channel single top-quark contribution is nor- 
malised to the observed cross section as measured from 
the combination of all four channels. 

6. Systematic uncertainties 

Systematic uncertainties on the normalisation of the 
individual backgrounds and on the signal acceptance af- 
fect the measured single top-quark f-channel cross sec- 
tion. In the NN analysis the shape of each individual 
prediction is also affected; both the rate and the shape 
uncertainties are taken into account by generating cor- 
related pseudo-experiments. The impact of the system- 
atic uncertainties on the f-channel cross-section mea- 
surement is estimated from these pseudo-experiments. 
The uncertainties can be split into the following cate- 
gories: 

Object modelling. Systematic uncertainties due to the 
residual differences between data and Monte-Carlo sim- 
ulation for the reconstruction and energy calibration of 
jets, electrons and muons are propagated in the analysis. 
The main source of object modelling uncertainty comes 
from the jet energy scale (JES), including the modelling 
of pile-up, as well as b-jet identification. Other com- 
ponents include lepton energy scale and lepton and jet 
identification efficiencies. The JES uncertainty has been 
evaluated using 2010 data I'sT]. Additional contribu- 
tions to this uncertainty due to the larger pile-up effects 
in 201 1 data are included and range from less than 1% 
to 5% as a function of jet pj and rj. For Zj-quark jets 
a JES uncertainty of 0.8% to 2.5%, depending on the 
jet pf, is added in quadrature to the JES uncertainty. 
Scale factors, determined from collision data (3^, are 



applied to correct the ^j-tagging performance in simu- 
lated events to match the data. Both b-jets and c-jets in 
simulation use the same ^-tagging scale factors with un- 
certainties that depend on the pj and rj of the jet. The 
uncertainties on the scale factors vary from 10% to 15% 
for fe-quarkjets and from 20% to 30% for c-quark jets. 
For light-quark jets the mis-tagging uncertainty ranges 
from 20% to 50% as a function of jet pj and rj. Other 
minor uncertainties are assigned to the reconstruction of 
£'™'^^ and to account for the impact of pile-up collisions 
on £'™'*^ Finally, a systematic uncertainty was also as- 
signed to account for temporary failures of parts of the 
LAr calorimeter readout during part of the data-taking 
period, which was not modelled in the MC samples. 
Monte-Carlo generators and PDFs. Systematic uncer- 
tainties arising from the modelling of the single top- 
quark signal and the tt background are taken into ac- 
count. The largest contributions come from the mod- 
elling of parton showers and hadronisation, estimated 
by interchanging the modelling between PYTHIA and 
HERWIG, and from the amount of initial-state and 
final-state radiation (ISR/FSR), estimated using dedi- 
cated AcerMC samples interfaced to PYTHIA where 
parameters controlling the ISR/FSR emission are var- 
ied in a rang^ consistent with those used in the Peru- 
gia Hard/Soft tune variations fell . The uncertainty due 



The default PYTHIA values of these parameters are : PARP(67) 
= 4.0 and PARP(64) = 1.0 for ISR and PARP(72) = 0.192 GeV, 
PARJ(82) = 1 .0 GeV for FSR. To decrease (increase) ISR, the param- 
eters PARP(67) and PARP(64) are set to 0.5 and 4.0 (6.0 and 0.25), re- 
spectively. To decrease (increase) FSR, the parameters PARP(72) and 
PARJ(82) are set to 0.096 GeV and 2.0 GeV (0.384 GeV and 0.5 GeV), 
respectively. Samples of simulated events are produced with six dif- 



8 




Lepton Charge and Jet Number m(lvb) [GeV] m(lvb) [GeV] 



(a) (b) (c) 



Figure 4: (a) Distribution of the lepton charge after the full cut-based selection for 2-jet and 3-jet events. (b,c) invariant mass of the fo-tagged jet, 
the charged lepton, and the neutrino, m(tvb), for the fo-tagged sample for 2-jet (b) and 3-jet (c) events after applying all cut-based selections except 
for the cut on m(fvh). In all three distributions the (-channel single top-quark contribution is normalised to the observed cross section obtained with 
the cut-based analysis. The last histogram bin includes overflows. 



to the choice of the single top-quark f-channel signal 
generator is estimated from the difference between Ac- 
erMC and MCFM predictions (4^. The modelUng un- 
certainty for the tt background is evaluated by compar- 
ing the generators MC@NLO and POWHEG Slii] 
(with HERWIG showering). For the W-i-jets back- 
ground a shape uncertainty is assigned based on the 
variation of the choices of the matching scale and of 
the functional form of the factorisation scale in ALP- 
GEN. Systematic uncertainties related to the parton dis- 
tribution functions are taken into account for the signal 
and for all background processes which are modelled 
by simulated events. In addition to the nominal PDF set 
the alternative MSTW2008nlo68cl [47] and CTEQ6.6 
PDF sets are investigated. Events are reweighted ac- 
cording to each of the PDF uncertainty eigenvectors and 
the total uncertainty is evaluated following the proce- 
dure described in Ref. Ii36] . An additional uncertainty 
is assigned for the mis-modelling of jets in the forward 
\t]\ regions. A weight function is derived from the pre- 
tag sample by dividing the observed I77I distribution in 
data by the distribution obtained from simulated events, 
for 2-jet and 3-jet events. The event weights defined in 
this way are then applied to all simulated samples in the 
^j-tagged data set. The systematic uncertainty is derived 
from the one-sided difference between the weighted and 
the nominal samples. The impact of using simulation 
samples of limited size is also taken into account. 
Theoretical cross section normalisation. The tt, single- 
top quark Wt- and .s-channel backgrounds are nor- 
malised to their theory predictions with theoretical un- 
certainties of 7% and 4%, respectively H [loi 



ferent sets of pai'ameters settings: ISR increased (decreased), FSR 
increased (decreased), and a simultaneous increase (decrease) of ISR 
and FSR. 



1 in . The uncertainty on the diboson background is 
5% (3?]. 

Background normalisation to data. The multijet back- 
ground estimate has an uncertainty of 50%. The NN 
analysis places an uncertainty of 50% on the rate of 
events with W-i-heavy flavour jets and 30% on the rate 
of W-Hlight jets events. These uncertainties are used 
as constraints on the predictions when simultaneously 
determining the W-i-jets rates and the signal cross sec- 
tion. The cut-based analysis does not apply a global un- 
certainty on the W-hheavy flavour and W-nlight flavour 
rates, but considers separately the impact of the domi- 
nant sources of uncertainty on the data-derived W-i-jets 
nomalisation factors. This treatment allows the correla- 
tion between each component of uncertainty on the nor- 
malisation factors and the uncertainties on the W-t-jets 
rates to be taken into account. The Z-i-jets background 
normalisation has an uncertainty of 60%. 

Luminosity. The uncertainty on the integrated luminos- 
ity is 3.7% ElliS. 

Table [3] shows the contribution of each source of un- 
certainty to the total uncertainty on the measured t- 
channel cross section (AcTobs/cTobs) for the neural net- 
work analysis and for the cut-based analysis. 

7. Cross section measurements 

Both the cut-based and neural network analyses em- 
ploy a maximum-likelihood fit method to measure the 
single top-quark f-channel cross section. The general 
likelihood function is given by the product of the Pois- 
son likelihoods in the individual channels. The back- 
ground rates are constrained by Gaussian priors. We 
use the following equations: 
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Table 3: Breakdown of the contribution of each source of uncertainty 
to the total uncertainty of the measured (-channel cross section in data 
for the NN analysis and the cut-based analysis. Theoretical uncertain- 
ties are included in the "Other backgrounds" uncertainty category. 





AcTobs/o-obs [%] 


Source 


NN 


Cut-based 


LJdld Statistics 


±j 


±0 


Object modelling 






Jets 


±0 


+ j/— 4 


Z^-ttigging efficiency 


+ 13 


±12 


Misttigging rate 


+ 1 


+1 


Lepton 


±2 


±4 


£'™■'*^ calorimeter readout 
1 ' 


±2 


±2 


Monte-Cai-lo 






PDF 


+3 


±4 


Generator 


+4 


±7 


Parton shower 


+5 


±11 


ISR/FSR 


+ 14 


+ 19/-18 


Forward jet modelling 


+6M 


+7/-5 


MC statistics 


±3 


±4 


Background normalisation 






Multijets 


+4 


±2 


Other backgrounds 


+ 1 


±6 


Luminosity 


+4 


±4 


Total systematic uncertainties 


+24/-23 


+30/-27 


Total uncertainty 


+24 


+31/-28 



k is given by cc^ and ajk for signal and background re- 
spectively. For the NN-based analysis the sets of and 
Ojk constitute the probability densities (shapes) of the 
NN discriminants. The scale factors for signal and yS^ 
for the backgrounds are the parameters of the likelihood 
function that are fitted to the data. The Gaussian con- 
traints on the background scale factors, Aj, are set to the 
theoretical cross section uncertainty for all background 
processes that do not use data-based estimates (ff, Wt, 
i-channel and diboson), while W+jets backgrounds are 
constrained within their data-derived uncertainties. The 
multijet background is fixed to the value estimated from 
data. 

The systematic uncertainties on the cross section 
measurement are determined using a frequentist method 
based on a large number of pseudo-experiments. For 
each pseudo-experiment the expectation values of the 
backgrounds vj and of the signal Vs as well as the rel- 
ative distribution of events across the channels (a^ and 
ajk) are varied including all sources of uncertainties de- 
scribed in Section |6] and the f-channel cross section is 
measured with the maximum-likelihood fit. The distri- 
bution of measured cross sections is an estimator of the 
probability density function of all possible outcomes of 
the measurement and it is used to estimate the uncer- 
tainty on the actual measurement. 

The NN-based analysis yields, from a simultaneous 
measurement in the 2-jet and 3-jet channels, a cross sec- 
tion of 



cr, = 83 ± 4 (Stat.) (syst.) pb 



83 ± 20 pb. 



k=l 



\\G{J3)-l,Aj) with 



and n\ ^/3)-vj-ajk. 



where M is the number of channels and B the number 
of background processes. The cut-based analysis uses 
M - 4 channels separated by lepton charge and the 
number of jets, while in the NN-based analysis M is 
equal to 28, namely the number of bins of the NN dis- 
criminant in the 2-jet channel plus the number of bins of 
the NN discriminant in the 3-jet channel. Here Vs and vj 
are, respectively, the predicted number of signal events 
and the number of events of background j in the selected 
data set. The number of observed (expected) events in 
channel k is denoted by rik (fik)- The expected number of 
events in channel k is //^ for the signal and /i^^ for each 
background /'. The fraction of events falling in channel 



The significance of the observed signal corresponds to 
7.2 standard deviations (6.0 expected). This is com- 
puted using as a test statistic the Q-valne, which is de- 
fined as the ratio of the value of the likelihood function 
maximized for the Standard Model signal cross section 
to the value of the likelihood function maximized for 
zero signal. 

The cut-based analysis measures, by combining four 
different channels (positive and negative lepton charge, 
with two and three jets) a cross section of cr, - 
92 ;^25 pb, in good agreement with the NN-based mea- 
surement. The separation of candidate events according 
to the lepton charge allows individual measurements of 
the top-quark and top-antiquark cross sections, yield- 
ing the results o-(f) = 59+[^ pb and o-(f) = 33+[^ pb, 
that can be compared to the theoretically predicted cross 
sections of 41.9+^-^ pb and 22.7+','^, pb, respectively H. 

To test the compatibility, the two measurements from 
the NN-based and cut-based analyses are combined 
using the Best Linear Unbiased Estimator (BLUE) 
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method The correlation coefficient of the two anal- 
yses is 75% and was determined with ensemble tests 
including all systematic uncertainties. Based on the en- 
semble tests the two results are found to be compatible 
within one standard deviation. However, the combined 
result and its uncertainty for the observed cross section 
measurement does not significantly differ from the re- 
sult obtained with the NN analysis alone. 

8. Vtb measurement 

Single top-quark production in the f-channel pro- 
ceeds via a W-t-b vertex and the measured cross section 
is proportional to \Vtb\^, where V,h is the relevant CKM 
matrix element. In the Standard Model \Vtb\ is close to 
one, but new physics contributions could alter its value 
significantly. 

The \V,b\ measurement is independent of assumptions 
about the number of quark generations or about the uni- 
tarity of the CKM matrix. The only assumptions re- 
quired are that \Vtb\ » \V,d\, \Vts\ and that the W-t-b inter- 
action is an SM-like left-handed weak coupling. There- 
fore, the tt background rate is unaffected by a varia- 
tion of \Vib\ since decays to a potential higher genera- 
tion are prohibited by kinematics. On the other hand, 
rates of single-top quark Wt and .s-channel backgrounds 
also scale with IVf^P, but their contributions are small in 
the signal region that drives the maximum-likelihood fit 
measurement. The resulting variation on the total top- 
quark background yield is less than its systematic un- 
certainty and thus considered negligible. 

The value of IVrbf' is extracted by dividing the ob- 
served single top-quark f-channel cross section, mea- 
sured using the NN method, by the SM expecta- 
tion |@]. The experimental and theoretical uncertain- 
ties are added in quadrature. The result obtained is 
\Vrb\ = 1.13+°;[4(exp.) ± 0.02(theo.) = LB+jJ;!^. Re- 
stricting the range of \V,b\ to the interval [0, 1], as re- 
quired by the SM, a lower limit on \Vtb\ is extracted: 
\V,b\ > 0.75 at the 95% confidence level. 

9. Conclusion 

In summary, we present a measurement of the cross 
section of single top-quark production in the f-channel 
with the ATLAS detector in pp collisions at V* = 
7 TeV. The measurement is based on a neural network 
discriminant separating signal events from background 
and yields across section of 83+20 pb. The correspond- 
ing coupling at the W-t-b vertex is \Vib\ = l.O^JJ jj and 
the 95% confidence level lower limit on the CKM ma- 
trix element \ V,b\ is 0.75. 
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